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Abstract. Since February 1997 the BL Lac object Mkn 501 is in a "high state" 
' °f 7 _ra y emission. The HEGRA collaboration has studied Mkn 501 with their air 

Q\ , Cerenkov telescopes on La Palma. Here we report on observations with the 5 m 2 

1 telescope (threshold s=s 1.5 TeV) operated in a stand alone mode. We observed 

a rapidly varying flux between 0.5 to 6 times of that from the Crab Nebula. On 
average a Mkn 501 flux of (2 + 1.3 — 0.5) x 10 -11 cm~ 2 s -1 has been determined. 
£^ , The spectrum extends at least up to 10 TeV with an integral power law coefficient 

■ of 1.8 ± 0.2 and seems to be steeper than in 1996. 

S : INTRODUCTION 

CO 

The AGN Mkn 501 is known to be a VHE 7-ray emitter [1,2] while it has not 
been observed by EGRET in the HE domain above 100 MeV. It is a variable 
^ source with a significant increase in intensity since its discovery by [1] in 1995. 

Since about February 1997 the source is in a 'high state' of emission showing 
strong variability [3]. The HEGRA collaboration has studied Mkn 501 since 
early March with their air Cerenkov telescopes (ACT). The ACTs are part of 
a detector complex for the study of cosmic rays over an energy range between 
~ 500 GeV and 10 17 eV. The installation is located on the Roque de los 
Muchachos, La Palma, Mas Canarias (28.8° N, 17.8° W, 2200 m a.s.l.). In 
total, 6 ACTs (5 of 8.5 m 2 and one with 5 m 2 mirror area) are in operation. 
Four of the 8.5 m 2 ACTs were operated in the so-called stereo-mode while the 
5 m 2 ACT (hereafter called CT1) was operated in a stand-alone mode. Here 
we report on observations with CT1 while the results from the stereo system 
are reported elsewhere at this symposium and first data are already submitted 
for publication [4]. Compared to the system, CT1 has a higher threshold (1.5 
vs. ~ 0.5 TeV) and only a 3.25° camera of 127 pixels, but was available for 



longer observation time starting at larger zenith angles. Also, the DAQ has 
a larger dynamic range for pulse height recording and observations can be 
carried out under the presence of moon light [5]. Technical details of CT1 are 
given in [6]. 



DATA SAMPLE AND ANALYSIS 

Since the 9th of March (MJD 50517) the source Mkn 501 has been observed 
continuously whenever weather conditions permitted it (~ 30 % of time lost). 
During periods of moonlight, observations were carried out with reduced pho- 
tomultiplier (PM) HV resulting in a 30-60 % higher threshold. Only during 
the days of full moon and the source being closer than 30° to the moon the 
measurements had to be stopped because PM anode currents exceeded safe 
limits. Observation times were typically 2-5 hours per night with raw-data 
trigger rates of 0.5 - 0.9 Hz (zenith angle dependent). The trigger condition 
was a coincidence of 2 pixels (of 127) exceeding each 17 photoelectrons (PE). 
After MJD 50568 the threshold was lowered to 13 PE by increasing the HV 
by 6 % resulting in a raw- data trigger rate of up to 1.8 Hz. 

The data were recorded entirely in the so-called tracking mode. The fil- 
tering and processing of the data and background determination followed the 
prescription given in [2]. After applying some filter cuts to remove accidental 
triggers etc., standard Hillas image parameters were calculated. For the se- 
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FIGURE 1. (a) The ALPHA distribution for on-source data taken during moonless night 
and background after all cuts except that on ALPHA. The angle ALPHA describes the 
orientation of the shower image and is expected to be small for 7-showers from the direction 
of the assumed source-position. At ALPHA < 10° there are 2554 excess events (40.4 cr). 
(b) like (a) but for the data taken at times when the moon was above the horizon. At 
ALPHA < 10° there are 455 excess events (12.4 <r). (c) Excess event distribution obtained 
by varying the assumed source position by up to ±1° along the RA and DEC axes (see 
text). 



lection of 7-ray candidates modified Supercuts, so-called dynamical supercuts 
[7] were applied. These cuts vary with energy (equivalent parameter: SIZE), 
shower impact parameter (equivalent parameter: DIST) and the zenith angle 
9. The cuts were at first optimised with Monte Carlo (MC) data and then 
fine-tuned on previously recorded data from the Crab Nebula. 

Fig. la and b show the ALPHA distributions for the total data with 9 < 30° 
which was used to produce the light curve shown in figure 3 together with the 
background determined from previously recorded off-source data. Fig. lc 
shows a "lego plot" of excess events versus assumed source position. The 
excess is centered within 0.02° (= tracking error of CT1) of the nominal Mkn 
501 position and has an angular spread of a = (0.10 ± 0.02)°. For a short 
period, CT1 was pointed 0.3° away from the source and the data showed the 
excess with the same shift and identical spread. The significance of the daily 
signal varied between 1 and 15 a depending on the flux and threshold while the 
excess rate varied between 2 and 60 h _1 . The background, however, remained 
stable (see fig. 2). 

RESULTS 

The rates were converted into a flux taking into account the collection area 
(determined by Monte Carlo simulations), change of threshold due to 9, in- 
crease of threshold due to the different HV settings and the dynamical super- 
cuts (loss of 7-ray events of ~ 40 %). Details will be given in a forthcoming 
paper. In order to assess the variability of the source, we present here the 
lightcurve for the data at 9 < 30° with an average threshold around 1.7 TeV 
(Fig. 2). The threshold rises with zenith angle and with reduced HV therefore 
the flux data from the concerned observations were extrapolated to 1.5 TeV 
by using a power law coefficient of -1.8 (see below). The complete analysis of 
the data up to 9 = 60° will be presented elsewhere. The daily flux showed 
large variations sometimes exceeding the CRAB flux by up to a factor 6. We 
determine an average flux (MJD 50517 to MJD 50608) of 

i ? Mkn501,March-Juncl99 7 (^ > 1.5TeV) = (2 + 1.3 - 0.5) X lO^CnrV 1 

where the errors given are mainly systematic. The comparable flux in 
1996 was [2] F M kn50i,M a rch-A ug usti996(£ > 1.5TeV) = (2.3(±0.4) stat (+1.5 - 
0.6) syst ) x lO-^cm-V 1 and the Crab Nebula flux [2] F Cra b,i996(£ > 1.5TeV) = 
(7.7(±1.0) stat (+4.6 - 1.9) syst ) x lO-^cm-V 1 . 

From MJD 50548 to 50554 (April 9 - 15, 1997) EGRET was pointed to 
the direction of Mkn 501 but no evidence for 7-ray emission above 100 MeV 
has been found from the preliminary analysis [8] with an upper limit in the 
range of the previous one of 1.7 x 10 _7 cm _2 s _1 while we measured during the 
same time a mean flux of (3 + 1.9 — 0.8) x 10 _11 cm _2 s _1 . This indicates that 
the spectrum must be very hard with an integral power law index between 



E 

o 



> 



A 



• nom. HV 
■ HV +6% 

O HV -4%, moon 

□ HV-8%, moon 

A nom. HV, moon 

* Sensitivity for 1 crh~ ,/2 

A Sensitivity for 1 ah" 1/2 , moon 



fit 



1 I 
I* 



HEGRA CT1 , Target! Mkn 50' 



1 f , 



1 . 



************ *** 



50520 



50540 



50560 



50580 



« 1 



50600 



50620 
MJD 



FIGURE 2. The preliminary daily flux measurements from March to June f997 from 
the data at zenith angle 8 < 30°. The different symbols indicate the different data taking 
conditions. Also indicated is the sensitivity, i.e. the flux corresponding to a 7-rate which is 
equal to the squareroot of the background rate in h _1 ; this gives a measure of the stability 
of the general observing conditions. 
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FIGURE 3. Inte- 
gral 7-candidate excess rates 
vs. energy threshold Eq for 
all data taken in 1996 (220 
h on-time) and 1997 (54 h 
on-time in the March-April 
dataset and 80 h in the 
May- June dataset) at zenith 
angles 6 < 51°. The bars 
are the statistical errors. The 
threshold was raised by ap- 
plying a zenith angle de- 
pendent cut on SIZE. The 
zenith angle distributions of 
the three datasets up to 51° 
are similar. 



100 MeV and 1.5 TeV < 1. In Fig. 3 we show the average integral 7-excess 
rates as a function of energy thresholds E for the first and the second half 
of the current observation period and for the data from 1996 (the 1996 data 
published in [2] were reanalysed using the dynamical supercuts which resulted 
in a 6.8 sigma signal with 494 excess events). These rates can be turned into 
fluxes by dividing by the appropriate expected rate which has to be determined 
from Monte Carlo studies. Preliminary investigations show that up to about 
10 TeV the 1997 spectrum can be described by a power law with an integral 
coefficient of 1.8 ± 0.2 which is within the errors compatible with the slope of 
1.58 ±0.51 of the 1996 data but the convergence of the rates at larger energies 
seems to indicate that the spectrum in 1997 is steeper than it was in 1996 and 
that the flux at energies above 10 TeV is essentially unchanged. This needs 
further investigations. Both the 1996 and 1997 spectra are incompatible with 
a large cosmological IR background calculated by [9] (and references therein) 
although a single source observation is not sufficient to rule out all possible 
scenarios. 

As Mkn 501 was also observed by the Whipple and CAT groups (see their 
contributions to this symposium) and both locations are at different longitude 
it is expected that soon nearly continuous data will be available covering a long 
period. The ongoing observations of Mkn 501 at energies from radio waves 
to pa 10 TeV will provide a unique opportunity to study the multi-wavelength 
behaviour of BL Lac objects. 
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